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Abstract 



We obtain the diagonal reflection matrices for a recently introduced family of dilute 
^ ■ lattice models in which the A3 model can be viewed as an Ising model in a magnetic field. 

, We calculate the surface free energy from the crossing-unitarity relation and thus directly 

(«— ^ ' obtain the critical magnetic surface exponent Sg for L odd and surface specific heat exponent 

for L even in each of the various regimes. For L = 3 in the appropriate regime we obtain 
the Ising exponent Sg = —y-, which is the first determination of this exponent without the 
use of scaling relations. 



There has been a recent growth of interest in the interaction-round-a-face (IRF) formulation 
C I of lattice models in statistical mechanics in the presence of a boundary ||l|, ^, |^, [s], ^ . IRF 

Q I models such as the restricted solid-on-solid (RSOS) Al models and the dilute RSOS Al 

models [^, ^ are particularly attractive as their solution in terms of elliptic functions correspond 
to off-critical extensions in which the elliptic nome p measures the deviation from the critical 
point p = 0. We refer to these models here as the Al models Q and the dilute Al models Q. 
In the Al models p is temperature- like while for the dilute RSOS models p is temperature- like 
for L even but is magnetic-like for L odd. In particular, the ^3 model, in the appropriate 
regime, can be viewed as a critical Ising model in a thermal field. On the other hand, the dilute 
A3 model, in the appropriate regime, can be viewed as a critical Ising model in a magnetic field 
[P. The singular part of the free energy of the A-^ model yields the bulk Ising specific heat 
exponent aj, = 0]. Consideration of the surface free energy yields the known Ising surface 



specific heat exponent ||l^, 11| = 1 |5|. On the other hand, the dilute A^ model provided a 



direct calculation of the bulk Ising magnetic exponent 6h = 15, first from the singular behaviour 
of the free energy and later from a calculation of the order parameters [12|. 

In this article we consider the dilute Al models with open boundaries and derive the surface 
free energy from which we obtain the magnetic Ising surface exponent (5^ = — ^ in the appro- 
priate regime of the A3 model. This is the first direct calculation of this quantity without the 
use of scaling relations. 

The dilute A^ lattice models |9| are RSOS models with L heights built on the A^ Dynkin 
diagram with a loop at each node. The non-zero face weights of the off-critical dilute Al models 
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satisfy the star-triangle relation [T^] and are given by 
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The crossing factors S{a) are defined by 
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and i9i(it), 'diiu) are standard elliptic theta functions of nome p 



t9i (tx, p) = 2p^/^ sin u n (l - V" cos 2n + (l - p^") 
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Four different critical branches are defined by B 
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This yields eight separate regimes, according to the sign of p. The magnetic Ising point occurs 
in regime 2 with A = 

The integrable boundary weights are represented by a triangular form with three spins 
[0, |2|, ^. For the dilute models, in accordance with the adjacency condition of the model, we 
can define 



K[ a 



unless |a — 6| = 0, 1 and \a — c\ = 0, 1. 
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These are to satisfy the boundary version of the star-triangle equation (the reflection equation) 
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ensuring the commutativity of the transfer matrix T(u) defined by the elements 
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where a = {oq, ai, • • • ,aj^} and b = {bo, 6i, • • • , b^} and 
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Note that the formulation of the refiection equation and transfer matrix given in is a special 
case of the general one given here. The above formulation does not incorporate crossing sym- 
metry of the bulk weights and is applicable to An \ n > 2 |l^ for which a crossing symmetry 
does not exist. 

Another important relation is the boundary crossing relation, which here reads 
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It can be easily seen that the reflection equation reduces to flve nontrivial equations with 
three distinct forms for diagonal boundary weights. These three distinct forms correspond to 
the generic equations in the higher rank Bn\An'^ cases, which have been solved in |15|. The 
diagonal solution we flnd is 
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where 



9a{u) =M^ + u)i?i(-A + e + u)M2aX + ^ + u)i9^{-2aX + ^ + u). 

Here the parameter takes the values ^ = —A/2 mod (^7r/2+m7rr/2) where ^ and m are integers 
and Qaiu) has been fixed by crossing symmetry, otherwise Qaiu) may be taken arbitrarily. The 
freedom in ^ due to the quasi-periodicity of the elliptic functions gives four distinct solutions. 
All of these satisfy crossing symmetry up to a gauge factor. As for the bulk weights, height 
reversal symmetry is broken by the boundary weights for L odd, ensuring that the nome p can 
again be regarded as a magnetic field. 

At criticality the known JC-matrix solutions for the A^2 ^ vertex model |]l^ are recovered via 
the usual passage from face weights to vertex weights. The integrable boundary weights of the 
dilute 0(n) loop model |jl^ are also recovered in this limit after transforming to the diagonal 
orientation ||5|. 

The fusion procedure has been applied to the dilute Al models, resulting in the construction 
of both su{2) Ijl^ and su{3) fused face weights from the face weights given in (|^). The 
su{2) fusion rule provides the functional relation llB] 



T(n)T(u + 3A) = /(n)+T(2)(u) 
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where T^'^\u) is the transfer matrix of the fused model with fusion level 2. For periodic 
boundary conditions the matrix function f{u) is given by 



f{u) = [p{u)ri 

where / is the identity and 
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The fusion procedure can be carried out in a similar manner for the open boundary system, 
as has been done already for a number of IRF models |jl], |5|, |6|. It is necessary to fuse the 
boundary face weights (|l^)-(|l3|). We find the matrix function f{u) to be diagonal with element 
(c, d) given by 



f{u\,d = u:-{u)L0+{u)p'''{u)/p{2u) 
where the boundaries contribute the factors w~(n) and a;j'(u), with 
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Height c (yd) is located on the right (left) boundary. 

We now turn to the calculation of the magnetic surface exponents 6s from the surface free en- 
ergy. The finite-size corrections to the transfer matrix T(n) are contained in T^'^\u). According 
to the spirit of ^, ^ the boundary crossing unitarity relation 



Tiu)T{u + 3X) = f{u)c,d 



(20) 
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for the eigenvalues is sufficient to determine both the bulk and surface free energies. The two 
contributions can be separated out by setting T{u) = Tb{u)Ts{u) and defining Tf, = k^^ and 
Tg = Kg- The free energies per site then follow as fb{u) = — logKfc(ii) and fs{u) = — logKs{u). 
Thus for the bulk contribution we have 

Kb{u)Kb{u + 3A) = p{u). (21) 



This relation has already been used to determine the bulk free energy M, 12 1 via the inversion 



relation method |13, 20 1. The critical behaviour as p ^ 0, obtained by use of the Poisson 




summation formula |13], is |§, 12 1 

r"" ^-"^-^ (22) 

L even ^ ' 

where the values of 5b and are listed for the different regimes in Table 1. 

To determine the surface exponents 6s and as it is sufficient to consider only the relevant 
contribution to the surface free energy, which does not involve the explicit form of the bound- 
ary weights. By making use of the boundary crossing relation and after taking a convenient 
normalisation, the appropriate contribution is obtained by solving 

T?i(5A - 2n)i?i(6A - 2u)t?i(5A + 2n)T?i(6A + 2u) 
T9?(5A)i?f(6A) 

We solve this relation under the same analyticity assumptions as for the bulk calculation ||8|, |l2| . 
Specifically, we set p = e~^, make the relevant conjugate modulus transformations, Laurent 
expand logKs(M) in powers of exp(— 27rii/e) and match coefficients in (pSj). In regimes 1 and 2 
we obtain 
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For regimes 3 and 4 the Ihs of relation (|23D needs to be modified to Ks{u)ks{u + 3A — tt) for the 
appropriate analyticity strip, with end result 
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Application of the Poisson summation formula yields 

/ pi+i/<5. L odd 
y L even 

as p — > 0, where the values of 5s and Os are listed for the various regimes alongside the bulk 
values in Table 1. As in the bulk case, these exponents are magnetic for L odd and thermal for 
L even. 

For L = 3 in regime 2 we obtain the value 5s = ~y- This result is in agreement with the 
prediction for the two-dimensional Ising model in a magnetic field using the scaling relations 



between bulk and surface exponents |21, 22, 23, 0]R The negative value indicates that the 



surface magnetisation diverges towards infinity as the magnetic field goes to zero. Thus a 
surface magnetisation in zero field does not exist. Similar behaviour is known to occur in the 



spherical model |11]. As expected, the Ising specific heat exponents ab = and = 1 are 



recovered in regimes 1 and 4 for L = 2. 



^Actually there is a misprint in the tabulated result for Ss in Binder's review 111 
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Table 1: Magnetic and thermal critical bulk and surface exponents of the dilute Al models, 
regime 12 3 4 



r 3L 
L-2 

2(2L-i; 
3L 



3(L+2) L-2 
L-2 3(L+2) 



-3(L+2) L+4 
L+4 L-2 



2(L+4) -2(L+4) 
3(L+2) L-2 



2(2L+5) -6 
3(L+2) L-2 



3r ^odd 

Lodd 
-2(L-2) r 

L+4 L even 
LTi ^ even 
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